Improved Semileptonic Form Factor Calculations in Lattice QCD 
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We investigate the computational efficiency of two stochastic based alternatives to the Sequential 
Propagator Method used in Lattice QCD calculations of heavy-light semileptonic form factors. In the 
first method, we replace the sequential propagator, which couples the calculation of two of the three 
propagators required for the calculation, with a stochastic propagator so that the calculations of all 
three propagators are independent. This method is more flexible than the Sequential Propagator 
Method but introduces stochastic noise. We study the noise to determine when this method becomes 
competitive with the Sequential Propagator Method, and find that for any practical calculation it 
is competitive with or superior to the Sequential Propagator Method. We also examine a second 
stochastic method, the so-called "one-end trick", concluding it is relatively inefficient in this context. 

The investigation is carried out on two gauge field ensembles, using the non-perturbatively im- 
proved Wilson- Sheikholeslami-Wohlert action with Nf = 2 mass-degenerate sea quarks. The two 
ensembles have similar lattice spacings but different sea quark masses. We use the first stochastic 
method to extract (D(a)-improved, matched lattice results for the semileptonic form factors on the 
ensemble with lighter sea quarks, extracting /+(0). 



I. INTRODUCTION 

Experimental measurements of heavy-light semilep- 
tonic decays, combined with theoretical input, can be 
used to extract the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements \V U b\, \V c b\, \V c d\, and \V CS \. The deter- 
mination of these matrix elements provides constraints on 
the CKM Unitarity Triangle and thus tests the Standard 
Model. Conversely, \V cc i\ and |V^ S | are known with high 
precision, and can be used to test the lattice techniques 
involved in calculating the decay rates of D mesons. Re- 
cent progress reports of current D meson semileptonic 
decay calculations are presented in Refs. [1-3]. The most 
recent study where systematic errors are taken into ac- 
count was presented in Ref. [4]. In this paper we study 
the efficiency of two stochastic propagator based alterna- 
tives to the Sequential Propagator Method used in all of 
the above calculations. 

The calculation we focus on is the semileptonic decay 
of a heavy- light pseudoscalar meson (H = D,D S ) to a 
light-light pseudoscalar meson (P = 7r, K). For processes 
H — > P the differential decay rate can be parameterized 
as a product of CKM matrix elements, perturbatively 
known quantities, and poorly known non-perturbative 
quantities, 



A 3/ V)|/ + (? 2 )| 2 , (1) 




CKM 



' form factor 



pcrturbativcly known 



where q 2 = (pn — pp) 2 is the squared difference between 
the initial and final state four-momenta and I € {d, s}. 
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The greatest source of uncertainty in the theoretical cal- 
culation of the decay rate is due to the non-perturbative 
interactions parameterized by the form factor f + (q 2 ). 

These interactions appear in the hadronic matrix el- 
ement {H(p H )\V^{q 2 )\P{pp)), where V M = ipd^i is a 
weak flavor-changing vector current and ip c is the charm 
heavy-quark and ipi is the d or s daughter-quark spinor. 
The matrix element can be parameterized as a linear 
combination of the vector, f + , and scalar, fo, form fac- 
tors, 

(H( PH )\V,(q 2 )\P( PP )) = 

{ P H+Pp-q{m 2 H -m 2 P )/q 2 } fi f + {q 2 ) (2) 

+ {q(m 2 H - m p)/q 2 }jo(q 2 )- 

On the lattice the matrix element is extracted from three- 
point functions of the following form, 



C 3 (T,r,;p ff ,q) = £e 



-'PH-x q-y 



x (0|Vvy 5 V> G (x, T)VvY M ^(y, t)^ 7 5^(0, 0)|0) 
_ £ e -PH-x e q. y Xr (AC 1 (0, 0; x, T)) 



x tsM-^x, T; y, t) 7M Mf 1 (y, t; 0, 0) 75 ) 



(3) 



where tp u is the spectator-antiquark spinor, M ~ 1 is the 
propagator for quark flavor q, and T is the time-slice 
at which the sink is fixed. The three-point function is 
shown schematically in Fig. 1. In the limit of large time 
separations Eq. (3) has the form, 



lim C 3 (T,t;p Hl q) 

T»t>0 



Zp Zh 
1Ep IE pi 



! H{p H )\V,\P{pp))e-^e 



-E P t -E H (T-t) 



(4) 
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FIG. 1: The three-point function, Eq. (3). Point-to-all 
propagators are used for the daughter-quark and 
spectator-antiquark. The heavy-quark propagator is 
constructed either by sequential or by stochastic 
all-to-all propagator techniques. 

The wave- function amplitudes, Zp and Zh, and energies 
can be removed by forming ratios with or performing 
simultaneous fits to the P and H meson two-point func- 
tions, 

C 2 P (t) ^ IgjV*** , C? (t) *»* lgtr*»« . (5) 

Currently, the approach of most lattice calculations of 
semileptonic decay rates of D mesons involves construct- 
ing the three-point functions using a so-called sequen- 
tial (or extended) propagator, the Sequential Propaga- 
tor Method [5]. As is shown in Eq. (3), this method 
uses two light-quark, point-to-all propagators M~ x and 
Mj~ , which may or may not be mass degenerate. It 
is routinely the case that these computationally expen- 
sive (expensive relative to the heavy-quark propagator) 
light-quark propagators have been previously generated 
to study other lattice quantities, and then stored. For 
example, the pseudoscalar pion mass is required to guide 
chiral extrapolations of many lattice quantities. 

To construct the desired three-point functions, the 
light-quark propagators are combined via the Sequen- 
tial Propagator Method with the heavy-quark propaga- 
tor generated "on-the-fly" . The point-to-all light-quark 
propagators connect the source, at (0,0), to the current 
at (t,y) and sink at (T, x). A point source is sufficient 
because of translational invariance. The heavy-quark 
propagator, M^^Tjy,*) = 75 [M c - 1 (y ) i; x, T)]t 75 in 
Eq. (3), must then connect all lattice sites at the vector 
current time-slice to all lattice sites at the sink time-slice. 
This is achieved by solving the heavy-quark action on the 
spectator-antiquark propagator sink at a given time-slice 
T, with a momentum and T matrix insertion. Although 
inverting the heavy-quark action with a k correspond- 
ing to the charm quark mass is inexpensive, relative to 
the calculation of the light-quark propagators, this ap- 
proach requires an additional 12 heavy-quark solves for 
every spectator-quark, D meson momentum, and/or D 



meson smearing. If several combinations of these param- 
eters are wanted, the cost of calculating the heavy-quark 
propagators becomes significant. In particular, at small 
lattice spacings the momenta are more finely grained, and 
therefore more three-point functions will have an appre- 
ciable signal, requiring many more heavy-quark solves if 
all data are to be utilized. 

Another approach is the One-end Method. This 
method combines sequential propagator techniques with 
stochastic propagators, through the so-called "one-end 
trick" as described for two-point functions in Refs. [6, 7]. 
The method has been used successfully in form factor 
calculations in light-light systems [8, 9]. This method 
provides an additional volume average and momentum 
projection which may greatly reduce the gauge noise in 
the correlators. However, additional propagators are re- 
quired for every source and sink momentum combination, 
suggesting it is computationally expensive relative to the 
Sequential Propagator Method. Nevertheless, there is 
the possibility that the three-point function noise is so 
greatly reduced that very few momentum points would 
be required in the q 2 — > extrapolation and a gain in the 
overall efficiency is realized. 

A more flexible method than the previous two that 
we also investigate in this paper uses an all-to-all [10], 
stochastically generated, propagator to connect the sink 
to the current. We will refer to this approach as the 
Stochastic Propagator Method. With this method all 
momenta and sink smearings are in principle available 
using only a single spin-color inversion. In addition, no 
new heavy-quark inversions would be required to calcu- 
late the three-point function using multiple (previously 
generated) spectator-quark propagators on the same con- 
figuration, as might be useful for chiral extrapolations in 
the valence light-quark mass. There is also the interesting 
possibility that correlators could be generated with mul- 
tiple current and sink time slices with no additional inver- 
sions. As we will show, the stochastic noise is significant 
and must be reduced by averaging over multiple noise 
vectors and partitioning [11-15]. In particular, time par- 
titioning is required for efficient reduction of noise, even 
when using other variance reduction techniques such as 
the hopping parameter acceleration (HPA) [16-19]. 

In the following, we study and minimize the stochastic 
noise appearing in the methods using variance reduction 
techniques. Comparing the errors and costs of the three 
methods in Sec. V, we determine the One-end Method to 
be inefficient but the Stochastic Propagator Method to be 
promising. We discuss when the Stochastic Propagator 
Method is more efficient than the Sequential Propagator 
Method, and finally apply the method to determine /+(0) 
on a single ensemble. This study is an extension of our 
preliminary results presented in Ref . [20] . 
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TABLE I: Simulation Parameters for the 2 QCDSF 
ensembles H and £ used. The scale was set from 
r = 0.467 fm, which was obtained by chirally 
extrapolating the nucleon mass times ro, TOjvro, to the 
physical point at j3 = 5.29 [21]. 



parameter 


ensemble heavy (H) 


ensemble light (£) 


P 


5.29 


5.29 


Volume 


16 3 x 32 


24 3 x 48 


^sea 


0.13500 


0.13632 


flight, valence 


0.13500 


0.13609 


^hcavy, valence 


0.11628 


0.12444 


a 


0.089 fm 


0.076 fm 




~ 930 MeV 


- 290 MeV 


valence 


~ 930 MeV 


~ 430 MeV 


m D 


- 2.43 GeV 


~ 1.93 GeV 



II. CALCULATIONAL DETAILS 

Two QCDSF ensembles with two degenerate sea quark 
flavors [22] were used in this study, with the parameters 
shown in Tab. I. The ensembles were generated using 
the Wilson plaquette action for the gluons and the non- 
perturbatively improved Shcikholcslami-Wohlert action 
for the fermions. Both actions have errors starting at 
0(a 2 ). We refer to the ensemble with relatively heavy 
sea-quarks as %, and to the one with relatively light sea- 
quarks as C. Wuppertal smearing [23] on top of APE 
smeared gauge links is performed on the interpolating 
fields in the two- and three-point functions, using smear- 
ing parameters that optimize the overlap with the light- 
light meson ground state. 

The ensemble with the heavier sea-quark mass, H, 
is significantly cheaper to work with, and so we per- 
formed some tests on it exclusively. The main results 
in this paper will, however, focus on the ensemble with 
the lighter sea-quark mass, C, since these parameters are 
much closer to the physical limit and the behavior of the 
noise should be influenced by the particle masses, lattice 
volume etc. As it turned out, the three-point functions 
generated on ensemble % had a noise behavior very sim- 
ilar to that on ensemble C 1 suggesting the influence of 
these simulation parameters on the noise is not large. 



III. SEQUENTIAL AND STOCHASTIC 
METHODS 

A. Sequential Propagator Method 

As stated in the introduction, the standard method for 
calculating D meson semilcptonic three-point functions is 
the Sequential Propagator Method. The sequential prop- 
agator provides a way to calculate a heavy-quark prop- 
agator that connects all spatial sites x at the sink time- 



slice T, to all sites y and t at the vector current insertion: 
the propagator M^ 1 that appears in Eq. (3). Algorith- 
mically this amounts to taking a single time-slice T of the 
spectator-quark propagator, M~ 1 (x, T; 0, 0), and pro- 
jecting it onto the sink momentum p# , the sink smearing 
Wgnk, and the sink T matrix. Finally, the heavy-quark 
Dirac operator is inverted on this "Sequential source" to 
get the sequential propagator, 

G(y,t;p H ,T;O,0) = 

^M- 1 (y,t;x,T)r^e J P"- x M- 1 (x,T;O,0) . (6) 
— v „ y 

X , ~ 

sequential source 

The sequential propagator, G(0, 0; p, T; y, t) = 
75[G(y, t; p, T; 0, 0)] t 75, can then be combined with 
the daughter light-quark propagator, M ; _1 (y, t; 0, 0), 
and appropriate gamma matrices to arrive at Eq. (3). 
This method requires 12 heavy-quark solves (one for each 
spin and color index) for each distinct sink momentum 
Ph, sink smearing TT sn k, sink time T, spectator-quark 
mass, and sink T matrix. The computational effort 
required for this procedure can become significant if 
many combinations of these parameters are needed, as 
would be required in partially-quenched chiral extrap- 
olations or when using the Variational Method [24, 25] 
for studying excited states. 



B. Stochastic Propagator Method 

With this method, the sequential propagator is re- 
placed by an all-to-all propagator, or more precisely 
with a time-slice-to-all (T-to-all) propagator 1 . We con- 
struct T-to-all propagators by generating noise vectors 
7?^(x, T), I — l,...,N, using complex Z 2 noise, with 
the property, 

^ E ^ (stof (*) - + ® (V Vlv) , (7) 

t 

where x = (x, T), z = (z,T), and i,j are combined spin 
and color indices. We then solve the linear system with 
the charm quark Dirac operator M c ^j(x,y) and source 
vector ?7^(x, T), for the solution vector ipf\y, t): 

E M c ,kj (z, T; y, t)$f (y, t) - n f (z, T) 

j,y,t , g . 

=> $f (y, o = E M c] k (y> *; z > T )^ ] {*, t) ■ 

k,z 



A true all-to-all propagator, where all time-slices have support, 
could be used. As we discuss in Sec. IV A, this turns out to be 
inefficient. 
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The average of the product of the source and solution vec- 
tors provides an unbiased estimate of the T-to-all heavy- 
quark propagator, 

1 J>? ] (y, t)nf ] (x, T) = M-^ (y, t- x, T) 
I " ' 

T-to-all 

k,z (9) 

V v ' 

Z 2 noise ccO(lfVN) 

where the stochastic error decreases with N, the number 
of noise vectors used. In the limit N — > oo, the exact 
T-to-all propagator is obtained. 

A stochastic estimate of Eq. (3) can be constructed 
by combining the spectator-quark and daughter-quark 
propagators and the T-to-all heavy-quark sources and 
solutions in the following manner, 

C 3 (T,i; P H,q) = 

- Jf E ■&(e- ipH - x r i M- 1 (O ) 0;x,T)r^M(x,T) 

xe t «V(y,*)rM I - 1 (y,t;O,0)) 
+ 0(1/VN) 

(10) 

where the appropriate propagator smearings must be ap- 
plied and = Tj = 75 and r = 7^. Although the error 
in the correlators generated using the Stochastic Propa- 
gator Method will be greater than that of the exact result 
from the Sequential Propagator Method, this difference 
is insignificant if the gauge noise dominates the stochastic 
noise. 



C. One-end Method 

An implementation of the One-end method for light- 
light three-point functions is introduced in Ref. [7], and 
can be directly adapted to heavy-light three-point func- 
tions. The approach begins with two T-to-all propaga- 
tors generated from the same noise vectors, and then 
uses the sequential propagator technique to complete the 
three-point function. We describe the case where the 
noise vectors, ry k , are used as the source at time-slice T 
for the heavy-quark and spectator-quark solves. 

The D meson momentum, p#, is first applied to the 
source vector and then a heavy-quark solve performed to 
obtain the solution vector, 

^i(y, t ; p H , T) = Mcjk iv, *! T)e^ ~r,f (x, T) . 

fe,X 

(11) 



A light-quark solve is then performed on the same source 
vector to get, 

V^(y',*') =J2M- 1 j , k ,(y',t';z,T)r ]k e }(z,T). (12) 

The pion momentum, pp, and 75 are then applied to the 
sink time-slice t' = of ipf) ul which is then used as the 
sequential source for the daughter-quark solves, resulting 
in the sequential propagator, 

Gf (y, t; pp, 0; z, T) = £ M^, (y, t; y', 0) 75 e Jpp - y ' 

j',fe',y',z 

xM-j, fc ,(y',0;z,T)^ 1 (z,T). 

(13) 

Forming the average of the product of Gj with ipj " 
(only the spin and color indices are transposed) and in- 
serting the appropriate current momentum q and current 
gamma structure 7sT results in an estimate for Eq. (3), 

C 3 {T,t;p H ,q) 

= e Jq ' ! 'e jpp - y 'e- J P«- x Tr(75M ( r 1 (x,T;y,i)r 
x,y,y' 

x M-\y, t; y', Q) l5 M^(y', 0; x, T)> 
^ E eiq ' y ^ lf (y. *; PH , ThsTGf (y, t- pp, 0; z, T) . 

(14) 

The potential advantage of this method is that the addi- 
tional volume average at the source (over y') and explicit 
projection onto p# could reduce the stochastic and gauge 
noise. 

There is another possible procedure for constructing 
the three-point function using the One-end Method. The 
T-to-all propagators could instead be used for both light- 
quarks and the sequential propagator technique for the 
heavy-quark. The Sequential and Stochastic Propagator 
Methods also have similar alternatives, where the sequen- 
tial or stochastic propagator can be used for the daughter 
light-quark and point-to-all propagators for the heavy- 
and spectator-quarks. However, with these methods it is 
obvious that the additional light-quark solves the alter- 
native procedure requires would make for an inefficient 
calculation. This is especially true in practice, where the 
point-to-all light-quark propagators commonly have been 
generated to calculate other quantities and stored. It is, 
however, unclear which procedure is most efficient for the 
One-end Method, and the required smearings, momenta, 
and quark masses in the final calculation influence this 
decision. It is also possible that one procedure results in 
significantly smaller errors than the other. When a dis- 
tinction is necessary we will refer to the procedure out- 
lined in this section as the One-end Sink Method, and the 
alternative procedure as the One-end Source Method. 
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IV. NOISE REDUCTION TECHNIQUES 

A. Partitioning and Hopping Parameter 
Acceleration 

The noise reduction methods we studied are partition- 
ing [11-14] and Hopping Parameter Acceleration (HPA) 
[16, 18]. With partitioning we decompose the noise vector 
sources into a sum of M subspaces in space-time, spin, 
and/or color, 



M 

711 — 1 



(15) 



so that the sources have support only on their respective 
subspaces. We then solve for the solution vectors ipm on 
each subspace. Our T-to-all propagator is formed by first 
summing over the M solution subspaces of a particular 
noise vector, then averaging over the product of the N 
source and solution vectors, 

1 £ V^(y, (x, T) « (y, t; x , T) . (16) 



Partitioning has the effect of explicitly setting many of 
the terms contributing to the stochastic error to zero. 
It has a computational overhead of M additional inver- 
sions per noise vector, but has been seen in some previous 
calculations to reduce the noise faster than ^= [11, 12]. 

Note that in the methods described thus far, 77^ had sup- 
port on only a single time-slice, corresponding to full time 
partitioning. This is not a strict requirement of the cal- 
culation: 7jM could have support on multiple or even all 
time-slices, allowing us to calculate the three-point func- 
tion with multiple sinks. We examine time partitioning 
in addition to other partitionings in the next section. 

HPA also seeks to set certain stochastic error contri- 
butions to zero, in particular the stochastic error terms 
close to the diagonal, which, in general, have larger am- 
plitudes. Writing the Wilson action M as, 



2kM = 1 - kD , 
and expanding the inverse in powers of k, 



(17) 



k-l 



M~ 



2nJ2( kI) y = 2k^(k£>) 1 + {KDfM- 1 , (18) 



i=0 



i=0 



we see that the first term on the right hand side should 
not contribute to the propagator for distances greater 
than k lattice sites, because D only connects nearest 
neighbor sites. These terms do however contribute to 
the noise when M _1 is stochastically estimated, due 
to the off-diagonal noise terms in Eq. (16). We can 
explicitly remove these contributions by using M~y = 
[(/c.D) fe M _1 ] , where \x — y\/a > k and a is our lattice 
spacing. 



For three-point function calculations HPA is poten- 
tially useful when multiple time-slices of the noise source 
vector have support. The stochastic error, in these cases, 
can be reduced by employing the HPA when the dis- 
tance in Euclidean time between the sink and current 
are greater than k, \T — t\/a > k. This has the effect 
of reducing the noise from the time-slices between T and 
t. Unfortunately, as we demonstrate in the next section, 
HPA does not reduce the error enough for multiple time- 
slice partitioning to be efficient in calculations of heavy- 
light semileptonic form factors . In the case of full time 
partitioning the HPA has no effect, because there are no 
time-slices between T and t that have support. 



B. Comparison of noise reduction techniques 

Due to the many additional light-quark solves required 
in the One-end Method, it is far more expensive to test 
noise reduction techniques using it rather than using the 
Stochastic Propagator Method. Therefore, we focus our 
efforts on minimizing the noise in the Stochastic Propa- 
gator Method, and make the assumption that the same 
techniques arc effective in the One-end Method. In the 
end we find the most efficient noise reduction technique, 
simple time and spin partitioning, provides a similar 
noise reduction in both the Stochastic Propagator and 
One-end Methods. 

We introduce the ratio of correlators we refer to as 
^o(?max)i which will be used as the starting point for 
comparing noise reduction techniques and the three dif- 
ferent three-point methods, 

C 3 (T,t;p = 0,q = 0) 



^(gLx) 



where lim V (q, 

T»i>0 



max 



cmc?{T-t) ' (19) 

• ^-(7r(0)|V |^(0)). We later 



form and study ratios of other correlators with different 
momenta and T insertions, but find Vo(<7 max ) to be par- 
ticularly convenient because it is the statistically cleanest 
correlator, with = p/j = 0. 

By comparing the relative errors of Vo((7 max ) con- 
structed with different combinations of time partitioning 
and the HPA we determined full time partitioning to be 
crucial for efficient noise reduction. For this test we use 
spin and even-odd partitioning in addition to time par- 
titioning, in order to reduce the number of noise vectors 
required to achieve a clear signal. The three different 
time partitionings we consider are: no time partitioning 
(with HPA and without HPA), partial time partitioning 
with HPA (where the noise vector source has support on 
time-slices T = 0,12,24,36), and full time partitioning 
(where the noise vector source has support on only the 
sink time slice T = 24). The results are presented in 
Fig. 2 for a single configuration from ensemble C, where 
the Vo((7 max ) presented has the pion source fixed at t = 
and the D meson sink at T — 24 for all partitionings. 

When comparing the different partitionings it is inter- 
esting to recognize that the correlators with no time and 



V :p = (0,0,0) (? = (0,0,0) 



• 


seo 


▲ 


seo-hpa 


■ 


seo-pt-hpa 


♦ 


seo-t 




FIG. 2: (color online) Vo^^ax) on ensemble £, using three different time partitionings, in addition to spin even-odd 
(seo). The pion source is at t = and the D meson sink at T = 24. The lower plot shows the relative error of these 
correlators. We show: no time partitioning (blue circles, seo), no time partitioning with HPA (black triangles, 
seo-hpa), partial time partitioning with HPA where the stochastic source vector has support on time slices 
t = 0, 12, 24, 36 (red squares, seo-pt-hpa), and full time partitioning (green diamonds, seo-t). Only the pion source is 

smeared. 



partial time partitioning can have multiple sinks T at 
no additional cost, — 24 and 0, 12, 24, 36 respectively, 
whereas the full time partitioned correlator has only a 
single sink T = 24. In an attempt to acknowledge the 
different statistical power available in each case we con- 
struct Vb(<7max) f° r each partitioning from different num- 
bers of stochastic vectors. For no time partitioning we 
use 96 vectors, partial time partitioning 16 vectors, and 
full time partitioning 4 vectors, the simplifying assump- 
tion being that each additional time-slice provides an er- 
ror reduction similar in magnitude to an additional noise 
vector. 

As demonstrated in Fig. 2, the noise of the correlators 
without full time partitioning is so much greater than the 
full time partitioned correlator (at fixed cost) that it is 
not necessary to quantify precisely the level of improve- 
ment due to having multiple sinks available: there is no 
way this can be competitive. In fact, the noise increases 
so quickly with distance from the stochastic source that 
further than time-slice £ = 15 the correlators without full 
time partitioning lose a clear signal. In the case of the full 



time partitioned correlator the noise is, however, nearly 
constant on all time-slices. Therefore we choose full time 
partitioning as the basis of our partitioning scheme, and 
note that this does not increase our cost relative to the 
Sequential Propagator Method, because that method is 
also limited to using a single sink time-slice to construct 
correlators. 

Starting from full time partitioning, we then tried all 
combinations of spatial (even/odd), color, and spin par- 
titioning, as shown in Fig. 3, again for Vb(<7„ lax ) on en- 
semble C using a single configuration with 100 noise vec- 
tors. The pion source is again fixed at t = and the 
D meson sink at T = 24 for all partitionings. The blue 
line represents the expected decrease of the noise with 
increasing the number of full time partitioned vectors. 
Perhaps surprisingly, none of the alternative partition- 
ing methods provide significant improvement over exclu- 
sively using full time dilution: spin partitioning on its 
own appears to have smaller errors, although the effect 
is small. This result is consistent over the other con- 
figurations we examined, causing us to choose full time 
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FIG. 3: The relative error, averaged over the time slices, of Vb(g„ lax ) for 7 different combinations of even/odd (eo), 
spin (s), and color (c) partitioning as a function of the number of solves N- lnvs required to calculate the partitioned 
correlator. All correlators are fully time (t) partitioned. These results are from a single C ensemble configuration, 

where only the source is smeared. 



p =(0,0,0) 9 = (0,0,0) 




FIG. 4: Total error at t = 8 for the fixed sink method, 
compared with the equivalent fixed current result using 
6, 12, and 24 stochastic vectors. Source and sink are 
smeared. 



and spin partitioning as our preferred variance reduction 
technique for the rest of the paper. 

The last alteration of the calculation we examine here, 



which could have an effect on the noise, is placing the 
stochastic source at the vector current instead of at the 
sink. This would result in a correlator with a fixed cur- 
rent and variable sink position instead of a fixed sink 
and variable current. The gauge noise, and presumably 
the stochastic noise, is different for each approach. For 
this test we considered 375 configurations from ensem- 
ble W. We present results using both approaches on a 
single time-slice for Vo(9max)i where the 7r source is at 
time-slice 0, the current is at time-slice t = 8, and the D 
meson sink is at time-slice T = 16. For this particular 
set of source, current, and sink time-slices, which is in 
the middle of the plateau as shown later in the paper, 
the two approaches are identical up to the effects of the 
stochastic noise. These results are shown in Fig. 4 us- 
ing 6, 12, and 24 stochastic vectors. It can be seen that 
the fixed sink correlator always has smaller errors. The 
gauge noise line shown is from the Sequential Propaga- 
tor Method. It is encouraging to see that even with only 
6 stochastic vectors the stochastic noise is only about 
~ 30% of the total noise. 
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V. COMPARISON OF METHODS 



A. Cost estimation 



To compare the efficiency of the three methods a way 
of measuring the overall cost is necessary. We choose to 
measure these costs C in terms of heavy-quark solves 
required, and note that with our simulation parame- 
ters a light-quark solve Ci requires approximately 30 
times the computational effort of a heavy-quark solve Ch- 
Ci s=s 30 Ch The Sequential and Stochastic Propagator 
Methods will always require the same number of light- 
quark solves as each other, however the One-end method 
will require different numbers of light-quark solves de- 
pending on stochastic source location, momenta, smear- 
ings etc. required to build the data set. 

We use N fully time and spin partitioned stochastic 
vectors, each of which requires 4 solves to generate. The 
Stochastic Propagator Method generates correlators with 
all sink smearings, VF S nk, and all momenta, p^ and po, 
with a fixed number of solves, AN . It requires no ad- 
ditional heavy-quark solves for correlators with different 
spectator quark masses, m u . The Sequential Propagator 
Method requires 12 heavy-quark solves for every differ- 
ent parameter at the sink or different spectator light- 
quark mass. Of course different source smearings, W src , 
spectator light-quark masses, m u , and daughter light- 
quark masses, mi, each require 12 additional light-quark 
solves to generate each point-to-all light-quark propaga- 
tor. The One-end Method requires new light-quark and 
heavy-quark solves for any changes in the parameters of 
the calculation. 

In Tab. II we summarize how the costs of the meth- 
ods break down according to the different parameters 
required for a particular calculation, assuming the light- 
quarks are non-degencratc. The integer factor a parame- 
ter contributes to the cost is simply labeled by the param- 
eter in parentheses. For example, if 2 different pu were 
required for a calculation then (pr>) = 2. Ci and Ch label 
the cost of a light-quark and heavy-quark solve respec- 
tively. The Wuppertal smearings at the source (sink) can 
be divided over the propagator sources (sinks) in what- 
ever manner is most efficient. The effects of requiring 
the light-quarks to be degenerate are: halving the num- 
ber of light-quark inversions required for the Sequential 
and Stochastic Propagator Methods, reducing solves re- 
quired in the One-end Source approach by AN(m u )Ci 7 
and dividing the last term in the One-end Sink approach 
by a factor of (mi). Note that the costs of the Sequen- 
tial and Stochastic Methods only differ by the last term 
which is proportional to Ch- 



B. One-end Method 

First we compare correlators generated by the One-end 
Sink Method to correlators generated from the Sequential 
and Stochastic Propagator Methods. For the One-end 
Method we use only a single, full time and spin parti- 
tioned, stochastic vector per configuration. Only a single 
stochastic vector is necessary because the properties of 
Eq. (7) are valid across different configurations, as well 
as within a single configuration. Of course more stochas- 
tic vectors could be used, but as Tab. II shows the cost 
grows very rapidly with additional stochastic samples. 
For the relatively cheap heavy-quark stochastic propa- 
gator required for the Stochastic Propagator Method we 
use N = 24 fully time and spin partitioned vectors. The 
exact costs of the methods are entirely dependent on the 
needed final data set, but here we demonstrate that the 
One-end Method is, in any "realistic" calculation, far less 
efficient than the other 2 methods. 

In general we expect the stochastic error of the One- 
end Sink and Source Methods to be different, and in fact 
found the One-end Sink Method to have smaller errors. 
According to Tab. II the Sink procedure is likely to be 
much more expensive than the Source procedure, but we 
defer discussion of the cost until the end of this subsec- 
tion. 

A comparison of the three methods on 250 configura- 
tions can be seen in Figs. 5a-5c on ensemble £, where all 
correlators are smeared at the source and sink. Fig. 5a 
shows Vo(g^ ax ) for the three methods. In this figure it 
can be seen that the Stochastic Propagator and Sequen- 
tial Propagator Vo(qf nax ) have smaller errors than the 
One-end Vo(q^ naxi ). I R Fig- 5b, a plot of the Vi(|p| = 
0, |q| = 1) correlator where a single momentum is used, 
shows the correlator for the One-end method to have the 
smallest errors. This might be expected due to the ex- 
plicit momentum projection at the source. This is too 
naive a comparison though, because the other two meth- 
ods have 5 more rotationally equivalent correlators avail- 
able at negligible cost, which can be averaged. Fig. 5c 
shows the effect of averaging the 6 correlators avail- 
able when using the Stochastic and Sequential Methods. 
The One-end generated data, which have no rotationally 
equivalent correlators to average over, now have slightly 
larger errors than the averaged data from the other two 
methods. 

We have performed similar comparisons for all other 
momentum combinations we are interested in and also for 
selected momentum combinations on ensemble "H, find- 
ing similar results in all cases. Considering there is no 
reduction in errors using the One-end Method and its 
high cost in a realistic calculation, it is obvious that the 
One-end Method is far less efficient than the other meth- 
ods. For example, consider a calculation with degenerate 
light-quarks where (m u ) — (mi) = 3, (pr>) = (p w ) = 3, 
(Wsrc) = (Wsnk) = 1, and d = 30C h : 
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TABLE II: A breakdown of the costs of each method. The integer factor a particular parameter contributes is 
labeled by the parameter itself, where the parameter labels are defined in the text. Each point-to-all and sequential 
propagator require 12 solves, and each of the N stochastic vectors require 4 solves. The approximate costs of the 
light-quark and heavy-quark solves are labeled by Q and Ch respectively. Ci ~ 30Ch with our quark masses and 

lattice volumes. 



method computational cost 

sequential 12(m;)(W S rc)Ci + 12(m u )C i + l2(m u )(p D )(W S nk)C h 

stochastic 12(m;)(W sf c)C; + 12(ra„)C( + 4NC h 

one-end (source) 4AT(p x )(mi)(W Br c)Ci + 4N(m u )Ci + 48N(m u )(p D )(W sn k)C h 

one-end (sink) 4N(p D )(W ank )C h + 4N(m u )Ci + 4SN{m u ){Pn){mi){W S rc)Ci 



Cost (sequential) = 12(3)(30 C h ) + 12(3) (3) (C/J =U88C h 

Cost (stochastic) = 12(3) (30 C h ) + 96(C h ) = 1176 C h 

Cost(one-end source) = 4(3)(3)(30 C h ) + 4(3) (30 C h ) + 48(3)(3)(C /l ) = 1872 C h 

Cost (one-end sink) = 4(3) (Ck) + 4(3)(30 C h ) + 48(3)(3)(30 C h ) = 133372C h . 



V : p = (0,0,0) g = (0,0,0) 




0.01 
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t 

(a) (color online) The upper plot shows the scaled temporal 
component of the vector matrix element, Vb, at q 2 = (^axi from 
the three different methods using 250 configurations. The lower 

plot shows the relative error of the three correlators. The 
Stochastic Propagator based three-point function is labeled by 
spin — 22 (blue circles) , the Sequential Propagator by seq (red 
squares), and the one-end by one — end (green triangles). 

FIG. 5: Error comparison of the three methods. 250 
configurations from Ensemble C are used. 



As can be seen in the calculation above we chose a very 
expensive One-end procedure to perform this calculation, 
the One-end Sink Method; however, even the One-end 
Source Method would be much costlier than the Sequen- 
tial or Stochastic Propagator Method. We mention again 
that the One-end Sink Method was chosen to minimize 
the error, and so the cheaper One-end Source Method 
would be even noisier. The One-end Methods at these 
quark masses and volumes are evidently not close to be- 
ing competitive with the other two approaches, and we 
drop these methods from further consideration. 



C. Sequential Propagator versus Stochastic 
Propagator 

In this section we compare the two more promising 
methods, as suggested by the plots and rough cost cal- 
culation from the last subsection. We perform a direct 
comparison of the Sequential and Stochastic Propagator 
Methods to determine when it is advantageous to use the 
Stochastic Propagator Method. 

As a first step in an effort to understand the behavior 
of the noise present in both methods we examine the ef- 
fects of smearing. Va{q^ sl&x ) with two different smearings 
is shown in Figs. 6a-6b for the two ensembles. We re- 
mind the reader that in this setting, in the limit TV — > oo, 
the Stochastic Method will yield exactly the same result 
as the Sequential Method. The results from 100 configu- 
rations on ensemble % are presented in Fig. 6a, and for 
220 configurations on ensemble £ in Fig. 6b. The be- 
havior of the gauge noise can be directly seen from the 
Sequential Propagator Method correlators. It decreases 
slightly with distance from the D meson sink regardless 
of smearing. The contribution of the stochastic error 
can be determined by comparing with this purely gauge 
error: for all smearings and both ensembles it is neg- 
ligible compared to the gauge noise near the stochastic 
source at the D meson sink, but grows mildly with dis- 
tance from the source. The effects of smearing on the 
gauge and stochastic noises are both mild, and smearing 
clearly helps to produce cleaner plateaus. 

We now compare the total errors of both methods using 
the five statistically cleanest combinations of correlators 
with finite momenta on ensemble C, where each combina- 
tion of correlators corresponds to a particular fo(q 2 ) and 
/+(g 2 ). When averaging over so many equivalent correla- 
tors it is more straightforward to extract the form factors 
from the appropriate combination of matrix elements and 



10 



VJ : p = (0,0,0) g = (l,0,0) 



V : p =(0,0,0) q = (0,0,0) 
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(c) average over Vi , i = 1,2,3 

FIG. 5: cont. Error comparison of the three methods. 
250 configurations from Ensemble C are used. 



kinematic factors, and then perform the average. The en- 
ergies needed for this construction are determined from 
bootstrap ensembles of two-point fits. These are then 
combined with the corresponding bootstrap ensembles of 
the three-point correlators. 

We assume only three D meson momenta are desired 
(on these lattices, D mesons with momentum |pu| > V2 
have significantly distorted dispersion relations): 



Pd = (0,0,0), (1,0,0), (1,1,0). 



(20) 



We average over all rotationally equivalent correlators to 
improve the statistics; however, there are many correla- 
tors with rotationally equivalent sink momenta which are 
not available to the Sequential Propagator Method but 
are essentially free for the Stochastic Propagator Method. 
These additional momenta are utilized in the stochastic 
case. We also consider the same momentum range for 
the pion (on these lattices and valence quark masses cor- 
relators with |p^| > V2 are too noisy to produce a clear 
signal) . Both methods have the same rotationally equiv- 
alent source momenta available. Figs. 7a-7e are repre- 
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(a) (color online) Vb(9m ax ) from 100 configurations on ensemble 

H. The Stochastic Method (purple circles) and Sequential 
Method (green squares) correlators on the left are not smeared. 
The Stochastic Method (blue circles) and Sequential Method (red 
squares) correlators on the right are smeared at source and sink. 
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(b) (color online) Vo(g^ al:ix ) from 220 configurations on ensemble 

C. The Stochastic Method (purple circles) and Sequential 
Method (green squares) correlators on the left are smeared at the 

source but not the sink. The Stochastic Method (blue circles) 
and Sequential Method (red squares) correlators on the right are 
smeared at source and sink. 

FIG. 6: We compare the effects of smearing on the 
Vo(9max) correlator on the two ensembles. The smearing 
dependence of the wave function amplitudes, Zp and 
Zh in Eq. (19), causes the difference in scale for the 
different smearings. 



sentative of the range of behavior of the relative errors in 
our data set. The number of correlators averaged over is 
shown below each figure for both methods. 

The Sequential Propagator Method uses 12 x 3 = 36 
heavy-quark inversions to generate the correlators. The 
stochastic propagator method uses (TV = 24) x 4 = 96. 
As suggested in Figs. 4 and 6a-6b, less vectors could be 
used without a significant increase of the error. However, 
as seen in the cost estimation in the previous section of 
a realistic calculation, choosing N — 24 is still cheaper 
(although modestly) than the Sequential Method. If we 
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FIG. 7: Comparison of errors from correlators 
generated using the Stochastic and Sequential Methods. 
The data in the above plots are generated using 668 
configurations. The spatial momentum corresponding 
to the q 2 of the data is shown in each sub-caption. The 
number of rotationally equivalent correlators available 
(and used) for averaging is listed after each method's 
name in the subcaptions. 



performed a naive fixed cost comparison with no assump- 
tions about what constitutes a realistic data set, only 
9 spin diluted vectors (requiring 36 heavy-quark solves) 
could be used, and the stochastic contribution to the to- 
tal error in the Stochastic Propagator Method would in- 
crease approximately by a factor of 1.6. As the stochastic 
error is a small percentage of the total error (the worst 
case appears in Fig. 6b in time slices 1-10, where the 
stochastic error accounts for ~ 5% of the total error) the 
Stochastic Propagator Method would still be competitive 
with the Sequential Propagator Method for most data 
points. 

We have examined all correlators with an apprecia- 
ble signal and draw the conclusion that with the mod- 
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FIG. 7: cont. Comparison of errors from correlators 
generated using the Stochastic and Sequential Methods. 
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(e) Stochastic - 6, Sequential - 6 p = (0, 0, 0), q = (1, 0, 0) 

FIG. 7: cont. Comparison of errors from correlators 
generated using the Stochastic and Sequential Methods. 
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est number of 24 full time and spin partitioned stochas- 
tic estimates the error in the stochastic correlators is at 
worst comparable to the error in the sequential corre- 
lators, while at best the additional data available re- 
sult in smaller total errors for the Stochastic Propagator 
Method. It is also clear that the Stochastic Propagator 
Method is much more flexible, because additional sink 
smearings, sink interpolating fields, and spectator quark 
masses do not add significantly to the cost of the calcu- 
lation. In the next section we use the Stochastic Propa- 
gator Method generated correlators to perform ag 2 ^0 
extrapolation. 



VI. MATCHING AND RESULTS 

In order to connect our results to observables of phe- 
nomenological interest we perform the matching and 
0(a) improvement of the vector current. The matching 
calculation takes the form, 



V™\q 2 ) = Z V [V^+aicvdvT^q 2 )] 



(21) 



where T M „ = V> c cr A1 „'0i is the tensor current with a^ v — 
|bv>7y]- The matching factor Zy is known non- 
perturbatively [26, 27], while the coefficient of the im- 
provement term, cy, is known to one- loop in perturba- 
tion theory [28]. We have examined the effects of 0(a) 
improvement and they do not change the behavior of the 
noise. 

It is interesting to note that a similar relation to that 
used in Ref. [2], where /o(? 2 ) can be directly related to 
the scalar current, holds to the order at which we are 
working. The observation the authors make in Ref. [2] is 
that the continuum Ward Identity, 



q^n\V^\D) = (m c - m ; )HS cont |A) , 



(22) 



is valid on the lattice up to a multiplicative matching 
factor. For the HISQ action [29], which they use, they 
obtain, 



Z v q»(n\V™\D) = (m c - m ; )<7r|S lat | £),) 



(23) 



This relation holds because the right hand side is renor- 
malization group invariant when the same action is used 
for both valence quarks. 

Due to the anti-symmetry of the tensor operator this 
relationship is also valid when using Wilson-type quarks, 
and allows us to calculate fo(q 2 )- 



fo(q 2 ) 



m c — mi 
m 2 D — m\ 



(S [at ) + 0{a 2 ) . 



(24) 



where /o(0) = /+(0). We extract /o(0) using this ap- 
proach as well as using the conventional approach with 
the vector current. 

As in the comparisons shown in Fig. 7, to extract the 
form factors wc first combine the appropriate three-point 



correlators, which are proportional to the vector matrix 
element, using the same bootstrap samples for all ro- 
tationally equivalent components, so that the resulting 
correlators in time are proportional to the form factors: 



C 3 {T,t;p H ,q) 



T>t»0 Zp Zh 



1Ep *2Eh 



/o,+e" 



-E P t -E H (T-t) 



(25) 

We then form the ratio from the bootstrap ensembles of 
our two-point and three-point functions, 



/o, 



Z-kZd 



C 3 (T,t; p g ,q) 
C„(t)C D (T-t) 



(26) 



and remove Z n and Zp, using our two-point fit results. 
It can be seen in Fig. 7 that the plateaus in the first 
two correlators are not clear. The lack of clear plateaus 
requires us to include excited states in our fits. To do 
this we determine the energy differences between the first 
excited state and ground state from two-state fits to the 
two-point functions. Using the energy differences we can 
perform a three parameter fit to extract the form factors, 



Po +Pie 



-AE w t 



+ P2e 



-AE D (T-t) 



(27) 



and obtain reasonable \ 2 values. Although plateau fits 
were suitable for the remaining correlators we found ad- 
ditional time-slices and wider fitting windows could be 
included using Eq. (27) for all correlators, and use this 
fit model for all form factor results presented here. 

Rcnormalized and 0(a) improved lattice results us- 
ing the Stochastic Sink Method for f+(q 2 ) and fo(q 2 ) 
calculated on the 24 3 x 48 ensemble are presented in 
Fig. 8. The form factors are plotted versus (roq) 2 , 
where r^ 1 « 0.422 GeV [21]. 530 configurations with 
24 stochastic vectors and 4 different time sources on ev- 
ery lattice were used. We also present results, in Fig. 9, 
for the same calculation using the relation Eq. (24). This 
is advantageous because no matching coefficients are re- 
quired, removing a source of systematic uncertainty. The 
two approaches to constructing the fo correlators are in 
all cases comparable within statistical errors; however, 
we find the statistical errors of the scalar current based 
data to be slightly larger. The q 2 — > extrapolation is 
performed using the Becirevic and Kaidalov parameteri- 
zation [30], 



fo(q 2 ) = 



2 _ c(l - a) 



1 - qy/3 ' 



Mi 2 ) 



c(l — a) 



(28) 

with q = q/m* D . m* D is the mass of the heavy- light vec- 
tor meson and c, /?, a are fit parameters. Reasonable x 2 
were obtained using (V^> or (S), with /+(0) = 0.593(19) 
and /+(0) = 0.603(24), respectively. The errors are sta- 
tistical only. Of course to extract physical results from 
any of these methods we would need to perform a chiral 
and lattice spacing extrapolation. We note that these re- 
sults are in agreement with the value /+(0) = 0.64(3)(6), 
obtained in Ref. [4], the most recent lattice calculation 
where systematic errors are estimated. 
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VII. CONCLUSIONS 

We have studied two stochastic based methods, the 
One-end Method and Stochastic Propagator Method, for 
calculating three-point functions in heavy-light semilep- 
tonic decays. The efficiency of these methods was opti- 
mized by minimizing the statistical noise, while taking 
into account the cost incurred by the different set-ups 
and variance reduction techniques. 

We found that the noise was reduced when placing the 
stochastic sources at the sink rather than at the current 
insertion of the three-point functions. This enables us 
to vary the position of the current but requires that the 
distance between source and sink is hxed, very similar to 
the Sequential Propagator Method. 

In terms of noise reduction techniques, full time parti- 
tioning was found to be essential, with no or marginal im- 
provement provided by additional combinations of color, 
spin, and even-odd partitioning. 

The One-end Method naively appears promising be- 
cause of the additional volume average it induces, which 
also has the effect of projecting out the desired momen- 
tum at the source. This creates the possibility of reducing 
the statistical errors, but at a high computational cost. 
In our analysis this method clearly fails to outperform 
the other methods, suggesting it is a poor alternative in 
heavy-light three-point functions. 

The Stochastic Propagator Method is competitive with 
the Sequential Propagator Method, and in many cases 
superior, due to the additional momenta freely available. 
The Stochastic Propagator Method is also more versatile, 
since correlators with different sink smearings, sink inter- 
polating operators, and spectator quark masses can be 
computed at negligible cost. This flexibility would allow 
calculations with the Variational Method and/or partial 
quenching to be performed with greater efficiency. This 
should further improve with finer lattice spacing, where 
the momenta are more finely grained. 

We have used the Stochastic Propagator Method to 
calculate the semileptonic form factors at many differ- 
ent values of q 2 's, using both the vector current and 
scalar current. Our data are parameterized well by 
the Becirevic and Kaidalov expressions. The results 
for the q 2 -> extrapolation, / + (0) = 0.593(19) and 
/+(0) = 0.603(24), calculated from the vector and scalar 
current respectively, are comparable to previous deter- 
minations. In the future we will extend this calculation 
to additional quark masses and lattice spacings, taking 
advantage of the method's flexibility to improve plateaus 
and reduce statistical errors. 
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